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I In the context of Af = 4 SYM, we compute the finite 't Hooft coupling A correction 

to the non-linear second-order hydrodynamic coefficient A3 from a Kubo formula based 
on fully retarded three-point functions using AdS/CFT. Although A3 is known to vanish 
in the infinite 't Hooft coupling limit, we find that the finite A correction is non-zero. 
We also present a set of Kubo formulae for the non-linear coefficients Ai^.3, which is 
■ more convenient than the one that has appeared recently elsewhere. 

M. 

1 Introduction 

The growing body of work on connections between gravity and hydrodynamics in the context 
of AdS/CFT seems to point towards an emerging picture: dynamics of small fluctuations of 
certain black-brane backgrounds is captured by the relativistic hydrodynamic equations [1] . 
On the other hand, AdS/CFT technology has been employed in studying strongly correlated 
phenomena, in particular understanding the strongly coupled quark-gluon plasma (QGP) 
created at RHIC. It turns out that hydrodynamics with negligible shear viscosity to entropy 
density ratio accurately captures the essential features of the state. Hydrodynamics is a 
classical effective field theory organized as a gradient expansion. When spatial gradients 
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in units of the mean-free-path are large, higher terms in the gradient expansion, beyond 
first-order become increasingly important, for instance in early-time evolution of the QGP. 
On general grounds, the stress tensor of a fluid can be decomposed as 



where Teq is the local equilibrium piece. It is a function of the fluid four-velocity and 
energy density e which itself is constrained by the equation of state. The extra piece U^ u 
is the dissipative part. In the present work it is truncated at the second order in derivative 
expansion. We work in the Landau-Liftshitz frame u^H^ u = 0, where the separation ([1]) is 
meaningful. The fluid velocity is normalized as u^vl 1 = — 1. 

At second order, there is a plethora of possible operators which are consistent with 
symmetries [21 [3]. In a conformal field theory however, which will be our focus here, the 
number of independent couplings reduces down to five transport coefficients: rri, k, Ai, A2 
and A3, in addition to the shear viscosity 7/ which already appears at first-order in gradient 
expansion. The ordinary effective field theory logic suggests the following form for on 
a general curved background [4J 



where the angular brackets are defined as = \^ a /\ vl3 (A a p + Ap a ) - ±A^A a/3 A*/3 

and the shear tensor = 2V^?W and vorticity £V iV = ^A fia A l 'P(V a up — V^Uq,) are 
written correspondingly using the angular brackets. Note that the coefficients Ai 5 2,3 are 
non- linear and necessary in order for the second order hydrodynamics to be consistent [1]. 

There have been intensive efforts in evaluating these transport coefficients in various con- 
texts, ranging from strongly coupled N = 4 super Yang-Mills theory (SYM) with infinitely 
many colors and infinite 't Hooft coupling A pQ to perturbative QCD [5j. Interestingly, for 
Af=4 SYM at infinite 't Hooft coupling, the coefficient A3 vanishes pQ . For a generic theory 
at arbitrary coupling, A3 is non-vanishing. See for example [6] for evaluation of A3 at zero 
coupling in scalar theory. 

One convenient method to derive transport coefficients from a microscopic description 
is to use the so called Kubo formulae based on hydrodynamics. Recently, a set of Kubo 
formulae for the second order transport coefficients appeared in [6], which extends the 
previously known Kubo formulaedto further include the non-linear coefficients Ai 2,3. These 
Kubo formulae are expressed in terms of fully retarded three-point functions of the stress 
tensor. In this work, first we present our own set of Kubo formulae for the second order 
hydrodynamics. They also involve stress tensor three-point functions, however compared 
to [6] a different set of the stress tensor polarizations are chosen. Our choice proves more 
advantageous in a holographic setting. As previously mentioned, A3 vanishes for J\f = 4 
at A = 00. In this paper, using the Kubo formula we propose for A3, we work out the 
next-to-leading contribution to A3 at finite 't Hooft coupling. The coefficient A3 turns out 
to be non-zero. This new result is one of the central motivations behind writing this note. 
As exercises along the way, we calculate fully retarded three-point functions in AdS/CFT 
using the gravitational dual description of N=A SYM theory at infinite as well as finite 

3 Kubo formulae for m and k, involving stress tensor retarded two-point functions were already known 
from the work of Baier et al [3] . 
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coupling. More recently, proposal for computing three-point correlation functions appeared 
in [7], which seems to be equivalent to the prescription [8j. The results we find can be 
checked against the existing results in AdS/CFT using other methods [HE]. This is a check 
on the Kubo formulae we propose ([§]), ([9|) and (fTUj) . This also checks the prescription of [7] 
for computing causal retarded n-point functions in AdS/CFT. 

The organization of the paper is as follows. In section two we define the fully retarded 
three-point functions. Section three is devoted to explaining the setup. The Bulk three-point 
vertices and computation of the back-reactions on the metric, necessary for calculating the 
three-point functions, are discussed in section four. Then we move on to the finite 't Hooft 
coupling corrections and compute A3 at next-to-leading order. We conclude the presentation 
with a discussion. 



2 Causal 3-point functions 

We will only be working with AA=4 SYM in four-dimensions (which implies a five-dimensional 
bulk) to avoid lengthly dimension-dependent expressions among other reasons. In what fol- 
lows, Greek letters denote the boundary spacetime coordinates. These coordinates will also 
be referred to as (t, x, y, z) interchangeably. Uppercase Latin letters M, N, ■ ■ ■ are reserved 
for the five-dimensional bulk spacetime. 

There are various ways of defining fully retarded three-point functions; one is through 
functional derivatives of the real-time generating functional W with respect to the source 
terms. Alternatively, one could define them without introducing sources and merely by 
performing a path integral. In general, the two definitions above will differ by contact terms 
[9] . The definition which is more intimately connected to the gauge /gravity correspondence 
is indeed to use the effective action. In [HE], the closed-time-path formalism was employed 
to define the generating functional W 



e W[h!M] 



xe -f^zC E l'f 3 (z)] e -if^^d i yC[^(y)M} 



where h\ and hi are the metric fluctuations (i.e., source terms for the stress tensor) living on 
the upper and lower closed-time contour respectively. See [6] for more details on definitions. 
It is common to choose a new field basis (the so called ro-basis) by defining average sources 
and operators h r = x \ and T r = Tl t T2 as well as difference variables h a = h% — hi and 
T a = T\ — Ti. Variation of the effective action with respect to h a yields T r 

-2i dW 

(TO) = ( y (3) 

where one sets h a = and h r = h\ = hi = h (the classical background value) only after 
performing the functional derivative. More generally, the fully retarded n-point function is 
defined as 



(0, 



(-i) n ~ 1 (-2i) n d n W (J) 
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where c.t. stands for possible contact terms. Before proceeding, let us spell out our con- 
vention for Fourier transformation 



G^' PA (0, y,z) = J d 4 pd 4 Q GZa aP ' P \p, qV^^ (5) 

where p ■ x = —pot + p z z. The calligraphic letters say p are defined as p = p/(2irT), where 
p is a momentum and T is the temperature of the system. 

We are now in a position to present the new set of Kubo formulae, involving stress- 
tensor three-point functions for the second-order hydrodynamic coefficients. Consider hy- 
drodynamics in flat space in local equilibrium. By perturbing the system through small 
background metric disturbances = ij^ + h^, one can create shear flow a^ v and/or 
vorticity . Solving the hydrodynamics equations of motion, namely conservation law 
for the stress tensor, one can obtain a consistent fluid velocity background as a function 
of h^ u and determine the expansion of (Tr U ) in powers of from the constitutive rela- 
tions. Field theoretically on the other hand, one can expand the response of (Tr U ) to the 
perturbations in terms of the fully retarded n-point correlation functions 



<Tr(*)> = Gr a (x)-\jd*yG^(x,y)Kp(y) (6) 
+1 [ d 4 yd A zGZ a a ^ x (x,y,z)h a p(y)h pX (z) + 0(h 3 ). 



The following set of (flat boundary) metric fluctuations can be consistently turned on and 
results in convenient Kubo formulae: {h xy (t), h xz (t), h yz (t), h tx (z), ht y (z)}. Let us work out 
A3 as an example 

<r**> = -^h ty d 2 z h tx - ~htxd 2 z h ty + ?fd z h tx d z h ty + o(h 3 ), u » = (1,0,0,0) + o(h), (7) 

where only contributions from h tx / ty are kept. Recalling the definition ([6]), one obtains 



A 3 = -4lhnJ pz d qz G2t' tv (P,q), « = .Hm^G«>,g). (8) 



Similarly, one can extract the following relations 

r, = i lim q), 2777V - « = lim ff), (9) 

p,q-+0 Oqo P,q^0 dpQ 

X 1 =r,r^- li m -^—G*y aa xz > yz ( P ,q), 

p,q->0 OpoOq 

where po + qo ^ in the Kubo formula for ij. We also have 

A 2 = 2I77W - 4 lim J^ G ^\p, q). (10) 

3 The setup 

In this section we will be concerned with the stress-tensor correlation functions in the infinite 
't Hooft coupling limit, thus it is sufficient to restrict ourselves to the two-derivative sector of 
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IIB supergravity action involving only gravity and a negative cosmological constant Later 
in the text, we manage to compute next-to-leading finite 't Hooft coupling contribution 
to A3. As mentioned before, this is interesting since the leading contribution is zero. To 
accomplish this task, one is required to go beyond the two-derivative action (llip . More 
explanations will appear in the section ([3]). For the present section the action, including 
counter-terms is 

A " f d 5 x^(R + + 2 [ d 4 xV^h K - % / d?xy/=R[. (11) 



The background solution representing thermal equilibrium in the boundary theory is a 
Schwarzschild-AdS black hole with a planar horizon. The following coordinate system is 
convenient 

ds 2 = g^ N dx M dx N = ^^[-fdt 2 + dx 2 + dy 2 + dz 2 } + -^du 2 , (12) 

u Au z j 

where Hawking temperature of the hole is T, L is the AdS radius and / = 1 — u 2 . The 
black hole horizon is at u = 1 and the boundary of AdS sits at u = 0. The stress tensor in 
asymptotically AdSs space is given by 

N 2 q/v 2 T 2 

!*» = -^(^ - h^K) - W (^ - T an, (13) 

where (T^ u ) = lim u _>o ■\[—hT iW and /i^ is the induced metric on the boundary at 

u = 0. The extrinsic curvature of the boundary is defined as K^ v = —\{V^n u + V^ra^), 
where is the unit vector normal to the boundary pointing outward and K is the trace 
of K^y. The Einstein tensor G^ u = R^ u — \g^i/R is built out of the induced metric on the 
boundary h^. 

3.1 Recipe for three-point functions 

From now on we set the AdS radius L = 1. There are multiple but equivalent ways to 
write down the bulk cubic vertices (or equivalently the full three point amplitude). One 
approach is to work out the on-shell action for the cubic metric fluctuations. This is the 
route taken for instance in |10| . In this approach, one has to pay great attention to the 
boundary terms generated in the process. However, there exists an entirely equivalent 
approach. Although completely straightforward, the next few lines will outline this. The 
functional derivative V3 = $ H (b)^(b)\ H (b) evaluated on-shell directly yields the three-point 
correlation function, where we have suppressed the indices, is the boundary data and 
lim u _j.o &9p\ = SHp\- Varying the action once leads to 

SS bulk = J d 5 x^£ MN [g PQ ]5g MN + J d^xV^hT^Jg^^ = (14) 

d b x^-g£ MN [G PQ ]8g MN + {T^)[g pX ]8H$ , 

where £mn is the equation of motion and Qmn is an arbitrary background. Observe that if 
the equations of motion for the metric perturbations are solved up to second-order in fluctu- 
ations, namely, £mn[q < ,8g] = 0(5g 3 ), the first term in (fT4|) will not contribute to V3 and 



4 We will not be considering loops. The action (fTTjl will contain the full set of relevant bulk vertices. 
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only two functional derivatives of the last term will survive. Since the last term is already 
a boundary term, no bulk integration over the radial direction ever shows up. This will be 
the definition for the three-point functions, we work with in this paper. As for the external 
wavefuctions, the usual boundary-to-bulk propagators (and their complex conjugates) are 
all that is required. We work them out in the next section. As emphasized in the full 
Kruskal-Szekeres extension of the black-brane background (to introduce Schwinger-Keldysh 
doubler fields) is irrelevant. The fully retarded functions we are concerned with here are 
causal quantities, therefore, on physical grounds only one causally-connected quadrant of 
the fully extended AdS-Schwarzschild spacetime should suffice for computing them. 



3.2 Solving the equations of motion with metric back-reaction 

As alluded to previously, in order to compute the stress-tensor three-point functions, it is 
required to evaluate the second-order back-reaction on the linearized metric fluctuations. 
This is formulated as 

9»u = 9 ( S N +h MN = g$ N+ eh^ N +e 2 h$ m £ M N[h W ,gW} = 0, £ MN [h^,g^] = S[(h (1 ^ 2 } 

(15) 

where e is a small parameter indicating the order in perturbation theory and £MN[h^ l \g^\ 
is the Einstein' s equation for the linear gravitons on the black-brane background. The 
term S[{h^) 2 } is the non-linear source stemming from interactions of the linear metric 
perturbations. Both h^ 1 ' and are taken to be incoming at the horizon (if is the 
largest-time insertion in the three-point function). At the boundary /tW = and = 
conditions are imposed, where is the boundary data for We demand regularity 
for all fluctuations at the horizon, although it turns out below that for certain metric 
polarizations demanding regularity at the horizon must be replaced by a stronger condition 
dictated by the equation of motion. The following metric polarizations are switched on 

h tx (u, z) = ^I^H tx (u)e^ z , h ty (u, z) = ^H^H ty ( U )e^ z . (16) 
u u 



h xz (t,u) = ^^H xz {u)e-^\ h yz (t,u) = ^I^H yz (v} ( -^ 

2 



hxy(t^Z^11^ — H X y(t : Z,Ti), 

where the four-momentum dependence of H xy is selected such that momentum is conserved 
for the cubic vertex of under study. This set of metric fluctuations is convenient as it can be 
consistently extended into the bulk to satisfy the bulk equations of motion and at the same 
time fluctuations are decoupled from each other. Let us focus on the three-point correlator 
Graa Z ' xz (p,q)- This will be relevant to Kubo formulae ([9]). The ansatz for H xy (t, z,u) for 
this correlation function is 

H xy (t, z, u) = H^e-' 1 ^^ 1 . (17) 
The equation of motion for H xy (u) including the non-linear source term is 



H 'U U ) + ^f^ H U u ) + (P \ + /2 q ° )2 g^) = ^H xz {u)H yz (u) + H' xz (u)H yz (u). (18) 
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It is easy to work out the equation of motion for H xz (u) 



2 2 

" -±^H' xz {u) + ^H xz (u) = 0, (19) 



where H yz satisfies the same equation but with qo replaced by po- Given that we are 
interested in fully retarded three-point correlation functions with xy-polarization as the 
largest-time insertion, we only need to include second-order back-reaction for H xy (created 
by a non- linear source in H xz j yz perturbations). Back-reaction on H xz j yz will contribute to 
other n-point functions, which we are not of interested in the present work. 

It will be sufficient to know the wavefunction for various metric polarizations as an 
expansion in frequency and momentum up to second-order. Let us define 

/C(u,ro) = (l-u)-** {1 - ln(l + u) (20) 

+to 2 [\ ln 2 (l +„) + (!-! ln(2)) ln(l + u) - ^dilog(i + |) + ^ - \ ln 2 (2)]}, 

then 

ffW po + qo), H x 1 J(u) = H x b JlC( U ,q ), H${u) = H$IC(u,p ). (21) 



The set of solutions (|2ip solves the first-order equations in (|15p . As for the back-reacted 
metric at the second-order, one obtains 

B%) = -H x b ]H^(l-ur^p oqo [\lnHl + u) (22) 

+ (1 - \ ln(2)) ln(l + u) - ~dilog(~ + |) + ^ - \ ln 2 (2)]. 

We need the renormalized holographic stress tensor (|13|) . expanded to the second-order in 
metric fluctuations. After substituting the solutions into the stress tensor and taking two 
derivatives with respect to the boundary data, one ends up with the following expression 
for the three-point function Graaa' vz (p,q) 

G x r y a r vz (P,q) = — i-ipo + q^N'T* (23) 

+i(l - ln(2))(p 2 + q 2 Q +p q )N 2 T 2 + O(p 3 ,p qlp 2 q , gg). 

Now let us turn our attention to Graa X ' ty (p z , q z )- H xy (t, z, u) will have the following form 

H xy (t,z,u)=H xy (u)e i ^ + ^ z . (24) 
The equation of motion for H tx / ty fluctuation is 

1 O 2 

H 'tx/ty( u ) - - H 'tx/ty( u ) + u ( u 2_ 1 ) g W") = 0, (25) 

where £2 stands for the momentum of the corresponding mode according to the definitions 
in (|16p . Note that there is an important subtlety regarding boundary conditions here. The 
equation of motion (125j) dictates a stronger boundary condition at the horizon; rather than 
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just asking for H tx / ty to be regular, one needs to set H tx j ty = at the horizon to obtain a 
regular solution. We need the following wavefunctions 



iJ« = flg> [1 - (p z + q z f In (1 + u)} + 0((p z + q,) 4 ), (26) 



and for H tx/ty 



«./* - - « 2 ) - -(i - «W + oca 4 ), (27) 

where Q stands for the spatial momentum of the mode according to (|16p . The back-reaction 
piece reads 

= H^H^[u 2 + (p z + q z f(u-u 2 -ln(l + u)) (28) 
+ (2u 2 - 3u + 2 In (1 + u))p,q 2 + 0(p 3 , p,q 2 , p 2 q„ q 3). 

evaluating the xy-component of the stress-tensor from (|13p . we get the following expression 
for the fully retarded function 

N 2 T 2 7r 2 /V 2 T 4 
« te '*(p„ 5.) = -^(P* + ^ 2 ) g— ■ (29) 

The momentum independent terms in (j29|) and (|23p can be understood in the following way. 
The xy-component of the energy-momentum tensor gives 

T xy = (e + P)u x u v + Pg xy + O(V). (30) 

where the error term denotes the gradient expansion. In addition, given we only care 
about second-order metric fluctuations, the first term in (|30p in zero frequency and spatial 
momentum limit will not contribute. As for the term proportional to pressure P, note that 
we have g xy = —h xt h yt + 0(h 3 ) (and g xy = h xz h yz + 0(/i 3 )), therefore 

lim lim G xy J x > ty (p, q) = -P, lim lim G?™*{p, q) = P (31) 

Similarly we show the A2 Kubo relation (|10p correctly reproduces the already known results. 
The relevant ansatz for H xy (t, z,u) is 

H xy (t, z, u) = H xy (u)e- ipot+i ^ z . (32) 

The solution for H xy , including back-reaction from H± y and H xz is straightforward to com- 
pute 

H xy = H® (1 - u)-«* [1 - (^Po + q 2 z ) In (1 + u)] - H®H® (1 - u)^p q. In (1 + u). (33) 
The xy-component of the stress-tensor yields 

G xy f z ' ty (Po,q z ) = ^Poq z . (34) 
Using values of the second-order hydrodynamic coefficients calculated elsewhere [TJ [U] 

vriV 2 T 3 n 2 -In (2) n x 7? In (2) 

^ = ^— ' T - = ^T-> Al = 2^T' X2 = —^-> As = °' (35) 

it is now trivial to see that, the Kubo formulae ([HD, ((UJ) and (|10|) . using the expressions for 
the three-point functions we have computed in ([23]) . ([29]) and ([MP reproduce the second 
order coefficients (f35j) . This concludes the warmup exercises. 
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4 Finite 't Hooft coupling corrections 



In this section we turn our attention to A3. Our goal will be to compute next-to-leading 
contribution to A3 (finite 't Hooft coupling correction). We need to consider a corrected 
IIB supergravity action. One known correction term is at order 0(a 3 ) and is quartic in 
the Weyl tensor. In [12], corrections involving F5 field strength were computed and further 
shown not to contribute to the corrected near extremal black D3-brane solution. In another 
work [14J, authors show that there is no need to work with the full 10-D action, if one is 
only concerned with thermodynamics and/or hydrodynamic response. Therefore we use the 
following a corrected gravity action 

N 2 r * 

S = ^ / cPxy/=g(R+ 12 + i3B), (36) 

where 

arr = 95^(2®^ -®^). b^ = V w) (37) 
'3 

C is the Weyl tensor and 7 = ^-£(3). The planar AdS-Schwarzschild, in the presence of 
the higher derivative term is corrected. The modified background looks like 

ds 2 = Wt { _ fZtdt 2 + dx 2 + dy 2 + ^2) + z Ju 2 

where 

Z t = l- 157(5n 2 + 5u 4 - 3u 6 ) + 0(j 2 ), Z u = 1 + 157(5u 2 + 5u A - 19u 6 ) + C(7 2 ). (39) 

The Hawking temperature associated to the hole is now 

T = T (l + 15 7 ). (40) 

According to the rules of AdS/CFT, world-sheet corrections proportional to 7 are 
mapped into finite 't Hooft coupling A corrections using a /L 2 = A -1 / 2 . Below, we find 
the 7 correction term to the fully retarded three-point function Graa X ' ty ■ Using the Kubo 
formula (|8|), we find corrections to k as well as next-to- leading piece of A3. The latter turns 
out to be non-vanishing. The correction piece to k has already been calculated elsewhere 
in [13]. Consequently, checking our result against that of [13] is an indirect check on our A3 
computation. Calculations are considerably tedious. Up to first order in 7, the solution for 
is given by 

H xy = H^[l-(p z + q z ) 2 ln(l+u) (41) 
+(P, + q,) 2 7(f ^ " 65n 5 + ^n 4 - ^n 3 + ^u 2 + 175 In (1 + u) - 175ti)]. 
The solution for H tx / ty reads 

H tx = ^ ) [(l-n 2 )-157n 2 (l-n 2 )(-3-u 4 + 5n 2 + 5)-p 2 n(l-n) (42) 
-57p 2 n 2 (l - u)(9u 5 + 25n 4 + u 3 + u 2 + u + 1)], 

and similarly for H ty with p z replaced by q z , where boundary condition H tx / ty = at the 
horizon was imposed. Given these first-order solutions, one can find the second-order back- 
reaction on H xy . We need to derive the H xy equation of motion by varying the action (f36|) . 
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including the non- linear quadratic source terms coming from H tx u y . We do not write the 
resulting equation of motion down, as it is extremely long and non-illuminating. We present 
the wavefunctions, as they may be useful for future explorers of the topic. To spell out the 
solution more efficiently, let us define 

H<® = ^^[F 1 + 7 F 2 + (p 2 + q z ) 2 F 3 -p 2 q 2 F 4 + 7(-p 2 q,G 1 + (p 2 + q,) 2 G 2 )] (43) 
where 

F 1 =u 2 , F 2 = 15u 2 (u 2 - l)(3u 4 - 5u 2 - 5), F 3 = u - u 2 - In (1 + it), (44) 

F 4 = -2ti 2 + 3-u-21n(l + u), (45) 

d = -90n 8 -^n 7 + 283u 6 -140w 5 + 135u 4 -— u 3 + 205u 2 -400u + 4001n (1 + «), (46) 
T 3 

„ o 7 283 c c 135 a 250 q 205 9 „ . , . 

G 2 = -45u 8 - 8(V + — u 6 - 65u 5 + — u A - — u 3 + — u 2 - 175u + 175 In (1 + u) . (47) 

_ £j O Zi 

At this point we must comment on the stress tensor in the presence of the quartic Weyl 
tensor contribution. We find that the stress tensor (|13[) is not affected by the higher deriva- 
tive correction in (|36f) . Indeed, contributions from the background proportional to 7 fall off 
too fast to make non-zero contribution to the stress tensor. So even in the presence of the 
term proportional to 7 in (]36p . the expression for the stress tensor remains unchanged. A 
similar observation was made in |15| . 

Here is the three-point function, including the finite A correction 

rwtxin, \ 7r 2 iV 2 T 4 15vr 2 iV 2 T 4 N 2 T 2 . 2 2 , 5iV 2 T 2 . 2 2 e . ^. 2 . 
G xy > tx > ty ( Pz ,q z ) = l+—^(p 2 z +ql) —(pl+ql-SPzqzh+Otf) 

(48) 

From this, for vanishing p z and q z , we recover the famous result of [16J, i.e., the finite 't 
Hooft coupling correction to the free-energy density of M = 4 SYM at finite temperature. 
Using the Kubo formula (jHJ) we obtain 

N 2 T 2 5N 2 T 2 
K = — ^ ( 49 ) 

which agrees with what appeared in [13]. Having passed these checks, we can read off A3 
using the Kubo relation (|HJ) 

A 3 = _^! 7 = - ( ^) 2 C(3)A-f N 2 T 2 . (50) 

This result, along with the set of Kubo formulae ([8]), and (jlOh are main results of this 
note. Note that the minus sign is just our convention used in writing the term proportional 
to A 3 in ((2j). 



4.1 Ward identities 

In this section, we provide further check on our A 3 result (150p . The coordinate invariance 
of the A/"=4 generating functional implies Ward identities. One particular Ward identity 
presented in [6], involves both two and three-point functions and only one momentum 

[^GlT^Q) + <* ")] + {^GZ M {Q) + (a O 0)] + ( 7 ^5) = 2G£T /3 ' 75 (Q,0). (51) 
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Note that we could always arrange for one four-momentum to vanish in the correlator 
relevant to the A3 computation; one could pick q z = —p z and Po = Qo = 0. From (15 1 [) we 
find 

G t r f x (p) = G^y(p,0), (52) 

where the fact that G tx ' tx (p) = G ty ' ty (—p) in our setup was used. Note that after Euclidean 
continuation of (|52p . one can reinterpret the polarizations and write H 

Gf tx ( Pz ) = G x E M (p z ,-p z ). (53) 

The correlator G xy > tx ' ty (p z , —p z ) obviously contains less information than the full correlator 
with the two momenta independent but it is easy to check (using hydrodynamics) that it 
is sensitive to a combination of k and A3 

G% ' tX ' ty (P*, -P,) = -P + (* + \h)p 2 z + 0(pl). (54) 

The Ward identity ([53j) indicates that G t ^' tx (p z ) should know about a combination of A3 
and k according to (|54h . We evaluate the Euclidean two-point function in (|53l) and find 
that indeed the Ward identity (|53p holds 



„tx,tx, \ n xy,tx,ty, s n ± i\ ± 2 oojv 1 2 ion jv 1 2 
G]3 UM = G E KPz> -Pz) = 5 H 5— Pz~{ 5 Pz + 5 )7+0(7 )• 



tt 2 N 2 T 2 N 2 T 2 9 ,35/Y 2 T 2 9 15vr 2 iV 2 r 4 , 

'E ' 

(55) 

This is yet another non-trivial test on our result, equation (|50p . 



5 Discussion 



In this section, we discuss implications of our findings. We also comment on future direc- 
tions. It was known from [I], that A3 = for the infinitely coupled TV = 4 SYM plasma. 
Our results indicate that at finite but large values of A, the coefficient A3 is non-vanising. 
In [6], it was noted that for generic weakly coupled field theories, A3 is non-zero. Although 
the same computation for N = 4 has not yet been completed [17J, it is likely that the same 
is true about N = 4 SYM at weak or zero coupling (perturbative expansion for A3 starts at 
<7°, for example in free scalar theory with one field [6]). 

Our findings in this paper are consistent with a picture where A3 of J\f = 4 SYM plasma 
starts off non-zero (or zero [17]) at vanishing coupling and approaches zero asymptotically 
as A -3 / 2 with a slope given by our result (f50|) at large values of A. The possibility of 
A3 = at vanishing coupling would be quite interesting, since our analysis would predict 
a non-monotonic behavior for A3 as a function of the coupling constant A. This has to be 
contrasted with, for instance rj/s which depends monotonically on the 't Hooft coupling. 
This remains to be seen and surely calls for further investigations. It would also be nice to 
find the non-linear second-order coefficients for rotating fluids in AdS / CFT along the lines 
of PSJ. 
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